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Abstract 

Based on the d'Alembert-Lagrange-Poincare variational principle, 
we formulate general equations of motion for mechanical systems sub- 
ject to nonlinear nonholonomic constraints, that do not involve La- 
grangian undetermined multipliers. We write these equations in a 
canonical form called the Poincare-Hamilton equations, and study a 
version of corresponding Poincare-Cartan integral invariant which are 
derived by means of a type of asynchronous variation of the Poincare 
variables of the problem that involve the variation of the time. As a 
consequence, it is shown that the invariance of a certain line integral 
under the motion of a mechanical system of the type considered char- 
acterizes the Poincare-Hamilton equations as underlying equations of 
the motion. As a special case, an invariant analogous to Poincare linear 
integral invariant is obtained. 



1 Introduction 

The theory of integral invariants (relative or absolute) introduced by Poincare 
in his pioneering work [19] and extended by Cartan[Jj has found many appli- 
cations in celestial mechanics and various other branches of science including 
analytical mechanics [3[51Eg[inilI21IIllI51IIEll221l2Ell2!Z!, mathematical 
physics [11], contact systems [23], fluid mechanics [23], optics and lasers [21] 
etc. In particular, the relative integral invariants are important for oscilla- 
tory systems. 
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It is known [18] that the motion of many mechanical systems is restricted 
by nonholonomic constraints so that every position in space is not accessible 
and in general, the variation of the velocity and velocity of variation are not 
identical. As a consequence, there are two viewpoints concerning the com- 
mutativity of the operation d (actual variation) and 5 (virtual variation) 
of kinematical quantities. The first viewpoint is due to Holder, Volterra 
and Hamel [16] which presumes that d — 5 operation's commutitivity is uni- 
versally applicable for all velocities whereas the second viewpoint due to 
Amaldi, Appell, Levi-Civita and Suslov p2J [H] assumes they commute only 
for independent generalized velocities. 

So far as the author knows, in all investigations concerning nonholonomic 
dynamical systems, the researchers usually use the first viewpoint in con- 
junction with the method of Lagrange undetermined multipliers and employ 
the Lagrangian variables or quasi-coordinates [9] to study the integral in- 
variant for such systems. 

In [JJ, the author has studied the Poincare-Cartan integral invariant for 
linear nonholonomic systems in terms of Lagrangian coordinates and in [2] 
author has employed the Poincare formalism based on the Lie theory of con- 
tinuous groups of transformations to generalize these results for holonomic 
systems. The aim of the present investigation is twofold: (i)to extend the 
Poincare formalism to include the generalized variation involving the vari- 
ation of the time and (ii) to examine the theory, without appealing to the 
method of Lagrange undetermined multipliers, from the second viewpoint 
[13j for nonlinear nonholonomic dynamical systems. To achieve this we fur- 
nish a precise definition of the virtual displacement in terms of the newly 
introduced parameters [2] corresponding to the asynchronous variation and 
obtain a generalization of the Poincare-Cartan integral invariant for such a 
systems. 



2 Generalized Variation 

Consider the motion of a holonomic dynamical system whose configuration 
at any time t along the actual (real) trajectory is determined by the set of 
C 2 functions 



of the time t of at least class C 2 and define the actual displacement in time 
'(it' by the relation 



x p (t) 



p = l s 2,..- 



n 



(2.1) 




p = l,2,... 



n 



(2.2) 
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Any one of oo ra configurations infinitely close to the actual trajectory and 
compatible with the constraints along a neighboring curve is defined by 

x*(t) = x p (t) + 5x p (2.3) 

The quantities 5x p , arbitrary functions of the time t of class C 1 , are the 
variations of the variables x p (t) that are obtained by keeping the time t 
fixed. The difference 

x*(t) - x p {t) = 5x p (2.4) 

is called the synchronous (or simultaneous) virtual variation of the variable 
x p (t) and is designated as 5-variation. 

Turning to a more general variation process, we define the position of the 
system in the actual (real) motion by the variables x p (t) and determine in 
a varied motion an infinitely closed position, by the function x*(t + At), 
where At is an infinitesimal change in time and is a differential function of 
the time t. 

Taking into account only small quantities of first order, we get 

x*(t + At) = x* p (t) + x* p At = x p (t) + 5x p + ipAt, (2.5) 

where the dot over the letter denotes differentiation with respect to time t. 
Introducing the symbol Ax p , we write 

x*(t + At) - x p (t) = Ax p = 5x p + ipAt, (2.6) 

which serves to define the variation Ax p of the variable x p where the time 
varies and the position along the varied path is not simultaneous to the 
actual path. These variations are called asynchronous (or non-simultaneous) 
virtual variations of the variable x p and are indicated by the symbol A- 
variation. Further, the A-variation of an arbitrary function F(x p ,t) of the 
class C 2 in the domain of the variables x p and time t can be readily obtained 
as 

AF = SF + FAt. (2.7) 

Since At is an infinitely small and differentiable function of the time t, by 
virtue of the relations (2.2) to (2.5), it can be easily shown that 

d(5xp)/dt = 5x p ■ (d(Ax p )/dt = Ax p + x p At), (2.8) 

which leads to the fact that the velocity of the variation and the varia- 
tion of the velocity are the same (not the same) according to the process 
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of 8- variation (A-variation), respectively. It is worth mentioning that the 
relation (2.7) play a vital role when the investigations are carried out on the 
basis of variational principles and their interpretation varies according to 
the dynamical system. It is to be noted that the relations (2.8), in fact, rep- 
resent the commutative (non-commutative) behavior of the operation d/dt 
and 8 {d/dt and A), respectively. 

We assume that the relation d{At) = A{dt) always holds for the indepen- 
dent variable t whereas it can be shown that A8x p = 8Ax p holds for the 
dependent variables x p . 

In order to extend the existing theory of integral invariants we transform 
the preceding analysis to the Poincare formalism that is based on the theory 
of continuous groups of transformations [20] . 

Let the group variables x p be holonomically connected to the cartesian co- 
ordinate u's by the relations 

Up = u p {x p ,t) (p,q = l,2,...,n) (2.9) 

and the Poincare parameters r] p (oj p ) of the actual (virtual) displacement by 
the relations [8 J 

7] p = A pq x q + A p (u p = A pq 8x q ) (2.10) 



where A pq , A p are functions of the x's and the time t and the matrix A pq is 
non-singular. In this paper the following conventions will be observed: (1) 
the summation convention is employed throughout the work; (2) the indices 
will have the range of the values A,/i, v = 0, 1,2, ...,n; p, q, r, s=l,2,...,.n; 
i, j, k,l=l,2,...,m < n; a, (3, 7, <5 = (m + 1), ...,n. 
We introduce the following: 

Definition 2.1. The change dG of an arbitrary function G(x p ,t) , during 
the actual displacement dx p in the time dt of the system, is determined by 
the relations |2j 



dG = [X G + r]pXpG]dt (p = 1, 2, n) (2.11) 

where Xq and X p , characterizing the infinitesimal displacement, are the 
operators defined by 
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which form a transitive group of operators if we require that the commuta- 
tors 

(Xq,X p ) = XqX p — X p Xq, (Xp, X q ) = XpX q — X g Xp, 
satisfy the relations 

(X ,X p ) = C q Qp X q , (X p ,X q ) = C r pq X r (p,q,r = 1,2,..., n) (2.13) 

Here Cq p , C pq , which depend on the x's and the time t, are the structure 

constants corresponding to the operators Xq,X p . 
We now give the following: 

Definition 2.2. In a simultaneous virtual displacement Sx±, 5x2, Sx n of 
the system, the change 5G in an arbitrary function G(x p ,t) is determined 
by the formula [8]: 

5G = u p X p G. (p=l,2,...,n) (2.14) 

Here the quantities u p are the parameters corresponding to the synchronous 
variation 5G of the function and are called the Poincare synchronous virtual 
displacement parameters or simply the virtual displacement parameters. 
In order to include the asynchronous variation AG of the function G(x p ,t) 
during an infinitesimal time At, we substitute from (2.11) and (2.14) into 
the formula (2.7) to obtain 

AG = {At)X G + (u p + T] p At)X p G 

If the virtual displacements are characterized by the asynchronous variations 
Ax p , then we need to extend definition 2.2. For this purpose, we use the no- 
tation xq = t, t]q = 1 and we write the expression (2.12) for the infinitesimal 
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displacement operators in the compact form: 

X^ = ^(xi,...,x n )-^-; £° = 1, e p = (/i,i/ = 0,l,2,...,n) (2.15) 

Further, we set 

ujq = 0, Axo = At = Qq, xq = rjo = 1 (2-16) 

and analogous to the Poincare parameters io v [8j, introduce the new param- 
eters corresponding to asynchronous variations by the relation 

Qfj, = + rj^Qo (2-17) 



which was first given in [U [2] . 

We now present the following new definition: 

Definition 2.3. The variation AG of an arbitrary function G(x p , t), during 
a virtual displacement Axi, Ax2, Ax n asynchronous to the time At = £lo, 
is described by the formula 

AG = %X tl G (m = 0,1,...,«) (2-18) 

where Qq, Q\, Q n , assumed to be functions of class G 2 , are the parameters 
of the virtual displacement corresponding to the asynchronous variation. For 
future reference we shall call them the asynchronous virtual displacement 
parameters. 
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Moreover, in view of (2.16-17), the formula (2.11) can be put in a more 
compact form as 

dG = [^X^G]dt (2.19) 

Since the system is holonomic, by virtue of relations (2.11-13) and the 
rule 5d = d5, the synchronous variations Sr) p of the parameters of real dis- 
placement are given by [8] 

SVp = ^ + C& q + (2.20) 

and, by means of the formula (2.7), it follows that 

Ar] p = 5r] p + fj p Q (2.21) 



Thus , using relations (2.17) and (2.20), it can readily be shown that 

a Vp = n p - r) p n + cp rVq n r + cg^n, (2.22) 

which represents the asynchronous variations At] p of the parameters of real 

displacement in terms of the new parameters fi's of possible displacements. 

It is known [T] that the process of ^-variation and integration commute 
with each other for holonomic dynamical systems, but from the preceding 
analysis it follows that the non-commutativity of the operations 'A' and 
differentiation 'd' implies the non-commutativity of the operation A and 
integration. That is, 

A [\)dt^ [ 2 A()dt (2.23) 
Jti Jti 

in general, even if the system is holonomic. Precisely, we state the following: 
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Lemma 2.4. Let J be a functional defined by the integral 

J= [fdt, (2.24) 
J o 



where f is an arbitrary function of rj p ,x p and possibly the time t. Then the 
asynchronous variation A J of the functional J is given by 

AJ = I (Af + fn )dt. (2.25) 
Jo 

For the proof of the lemma, readers are referred to However, one 
can see that the non-commutativity of the A-operation and the integration 
is obvious due to the presence of the quantity ftio in the expression (2.25). 

In a special case when the variation is synchronous, we have At = Qq = 
0. This, together with (2.17) implies that A = 5 and (2.25) includes the 
well-known result of the commutativity of (5-operation and integration for 
a holonomic dynamical system. But it is to be noted carefully that for 
nonholonomic system these results vary in a different way. See [13| for a 
discussion of this case. 



3 The two viewpoints 

From the analysis discussed in the preceding section, it follows that the 
commutativity (or non-commutativity) of the operations d with 5 (or A) 
is expressed by the transpositional relations (2.20) and (2.22) respectively, 
which play a crucial role in the dynamics of the nonholonomic systems. In 
fact, the study of integral invariants for such systems is intimately connected 
with these relations, since they involve certain variational procedure. In the 
sequel we briefly discuss the viewpoints about the validity of these relations 
which allow us to obtain a generalization of the theory of integral invariants. 

Let us consider the dynamical system whose configuration is determined 
by the group variables x p and moves subject to the nonholonomic constraints 
expressed by the (n — m) independent equations 

fa(r} P ,x p ,t) = (p=l,2,...,n;a = m + l,...,n), (3.1) 



S 



where the functions f a are not necessarily linear in the rj p 's and that the 
virtual displacement parameters u p 's satisfy the Chetaev's relations |8j 



drip 



u p = 0. (3.2) 



We further assume that the constraint equations (3.1) are expressible in the 
form 

fa(r) p ,x p ,t) = r] a - <p a (rji,x p ,t) = 0, (3.3) 
(p = 1, 2, n; i = 1, 2, m < n;a = m + 1, n) 
and the a/s satisfy the relations [13] 



uj a = -r—Ui (3.4) 
drji 

where we designate T] a (r]i) and u) a {uji) as the dependent (independent) Poincare 
parameters of real and virtual displacements respectively. 

In view of (3.2) or (3.4), we observe that the presence of nonholonomic 
constraints (3.1) or (3.3) leads to the fact that the w^'s are not independent. 
This reveals that the parameters of virtual displacement are not determined 
uniquely. Thus, there exists some arbitrariness in the determination of the 
derivatives Consequently, in using relation (2.20), we can adopt either 
of the two viewpoints given in the following: 

i) According to Holder, Volterra and Hamel, (2.20) may be used for the 
5- variation of all the parameters r] p of real displacement whether the system 
is holonomic or not. This, in view of the formula (2.14), allows us to express 
the variations 5f a of f a for the constraints (3.1) and (3.3) in the form [13j : 



Sfa 

and 



XpU + {c r 0p + c r qp v q )^ -j t \jb 



top (3.5) 



Sfa = 5l] a - S(j) a = A?LOi (3.6) 

respectively, by using the conditions (3.2) and (3.4). Here the quantities Af 
are determined by the relations 

A? = Jt (H) " ~ + + + + 

(q = 1,2, ...,n;i, j = 1, 2, m < n;a, (3 = m + 1, n) 

(3.7) 
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ii) According to Amaldi, Levi-Civita and Suslov as discussed in |13j . the 
relations (2.20) hold for nonholonomic systems. They may be used only for 
the variation Srji of the independent real parameters rji together with the 
assumption that the S- variation of the f a vanishes. Precisely, 

Svi = ^ + Ci q u q + C l qr r] q u r ; 6f a = (3.8) 

for the constraints of the type (3.3). 

The relations 5f a = enable us to obtain the synchronous virtual 
variations of the dependent ry's. To distinguish them from the variation 
5f(r]p,Xp,t) of arbitrary function of all the rfs, x's and t, we employ the 
symbol 5* to denote the variation of / when computed in terms of the in- 
dependent Poincare parameters by means of (3.3). Thus, from (2.14) and 
(3.4), it follows that 

S*f a = 5*rj a - 5*(f> a = (Af)* Uj (3.9) 
Here (A?)* are given by 

(3.10) 



where 





- ^Oj 


, n k U VP 
+ L ^d m 


trot 


- 


+ ^0/3 


K k - 






K a 




+ C q/3 


and 












X* 


= *3 











drjj 

9<f>f3 
drjj 



(3.11) 
(3.12) 



(3.13) 

4 Equations of Motion without Lagrange Multi- 
pliers 

We start with the general equation of dynamics expressing the d'Alembert- 
Lagrange-Poincare principle given by [13] 

d ( dL\ q dL dL 1 

dt{Wj- C ° P Wq - ^Wr ~ XPL \ UP = (41) 
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where L(r) p ,x p ,t) is the Lagrangian function that describes the dynamical 
behavior and the uj p determine the virtual displacement of the system. We 
assume that the system moves subject to the nonlinear nonholonomic con- 
straints (3.3) and adopt the second viewpoint. Since the a>'s are not inde- 
pendent and satisfy the conditions (3.4), the general equation of dynamics 
(4.1), as it stands, does not yield the equation of motion. Therefore, we 
write (4.1)in the form 



d di d<t> p d dL _ ^ p ^ p 

dt d Vj + d Vj dt d m { v + w d Vj ' 1 « + ^ d Vj )Vq 



- {XjL + ^XpL^u^O (42) 
(q = 1, 2, ...,n; i,j = 1, 2, ...,m < n;a, (3 = m + 1, n) 

Where we have used the relations (3.4) after separating the sum over p from 
1 to n into the sums over j from 1 to m and over (5 from (m + 1) to n. 
After expressing the r/^'s in terms of r/j's we call the resulting function 
L*(rjj,x p ,t). Thus, 

L*(r]i,x p ,t) = L(r] j ,(f) a (xp,r] j ,t),Xp,t), (4.3) 

as a result of operating by X p and differentiating partially with respect to 
rjj, we have 

dL dL dL* dL d6 a , . 

X p L* = X p L + -j—Xp^ — = -—-—-pL 4.4 
dri a dr]j drjj dr] a drjj 

Substituting from (4.3-4) into (4.2) and performing some mathematical ma- 
nipulations the general equation of dynamics takes the form 

' d ( dL \ f (r k +r k d ^ \ + (r k + r k d(f>f3 )n \ dL * 



(4.5) 



- ( Xj + ^X (3 )L* A* (X, + ^x r M. 



dr]j dr] a [ dt drjj dr\ 3 

, (r a , r a d( t>P ^, (ra , r oc d< t>P , 9(f> a ( k k d(j)f3 



{Cgj + C qf3-^-)Vq 







From the independence of the w's and (3.10-12), the last result leads to 

IC* - (Kfe + K ^ - x;L ' - { ^'0' ~ °- (46) 



11 



which are the required equations of motion independent of the Lagrange 
undetermined multipliers for the nonlinear nonholonomic dynamical system. 
These m equations determine the values of r]j(j = 1,2, ...,m) which, by 
virtue of the constraint equations (3.3), allow us to determine the values of 
the remaining (n — m) r/ a 's as functions of the x p 's and t. The values of x p 
are then calculated from the equations 

X p = X X p + 1] q X q X p (4.7) 

which are obtained from (2.11) by setting G — Xp. 

In order to obtain the canonical form of equations (4.6), we introduce 
the generalized momenta y p by the relations 

and we assume the transformations (4.8) to be invertible so that the ry p 's 
can be expressed in the form 

rjp = r)p(x q ,y q ,t) (p, q = 1, 2, n) (4.9) 

We note that (4.8) and (4.9) yield, respectively, y's and 77's as a linear func- 
tion of 77's and y's. 

Since all the 77's are not independent, it follows that y's are not indepen- 
dent as well. Therefore, we cannot work in the phase-space of 2n variables 
defined in [2], that is, 

Definition 4.1. A set of 2n independent quantities (x p ,y p ) is said to form 
the phase space if the conditions 

<57j = 0, u p ^0, 57]p^0, 5y p ^0 (p = l,2,..,n) 

hold throughout the motion of the system. 

To overcome this problem we introduce the new independent variables 
yj that are defined by means of equations of transformation 



dL* 
dm 



y*j = — ( 4 - 10 ) 
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which are assumed to be invertible and allow to express the m's as functions 
of y* , x p and time t in the form 

m = rjj(x p ,y*,t) (j = 1, 2, m < n;p = 1, 2, n) (4.11) 

Note that the constraint equations (3.3) are expressed in terms of the inde- 
pendent y*'s by means of the equations (4.11). 
We now give the following: 

Definition 4.2. The 2m quantities (xi, x m ; y*, y£J are said to form 
a 2m-dimensional reduced phase space provided that they are independent 
and satisfy the conditions 

St = 0; u)j ± 0; 5m ^ 0; Sy* + (4.12) 

throughout the motion of the dynamical system under the nonlinear non- 

holonomic constraints (3.3) with (3.4). 

Let us consider the motion of the system in this 2m-dimensional re- 
duced phase space and introduce the Hamiltonian function H*(x p ,y*,t) 
corresponding to Lagrangian L* 

H*(x p ,y*,t) = my* - L*(x p ,m,t) (4.13) 

Performing the 5- variation according to (2.14), in conjunction with (4.10), 
we have 

5H*{x p ,y*,t) = S[my* - L*{x p , Vj ,t)} (4.14) 
On one hand, using the second viewpoint and (4.10), we have 

5*H*(x p ,y*,t) = mS*y* - uj p X p L* 

Breaking the sum over the index p from 1 to n into the sums over the indices 
j from 1 to m and a from (m + 1) to n and using conditions (3.4) and (3.13), 
the last result assumes the form 

5*H*(x p , y*,t) = mS*y* - ujX*L* 
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which, in view of equation (4.6), becomes 

f) T * f)T 

6[r lj y*-L\x p , Vj ,t)}= Vj 6y*-u j [y*-(KK j +K k qjVq )^ ' 



On the other hand, we can find 



dr]k 3 dr] a 

(4.15) 



dH* 

5H*(x p ,y*,t) = -Q^-Syj +uj p X p H* 



which is equivalent to 



f) TT* 

5H*(x p , y*,t) = —5y* + UjX;H* (4.16) 

where we have used (3.4) and (3.12). 
From (4.15) and (4.16), it follows that 

f)H* f)T* r)T 

(— - rity + [X*H* + y* - (K^ + K k qjVq )— - {A«T{—T]u 3 = 

since all the 5y* and u>j are independent, the coefficient of each of the 8y*j 
and ojj vanishes and we obtain 

r) J-f* f)T * r)T 

m-Sj? « = -^ + (^ + ^*)^-^T(^)* (4l7) 

(j,k = 1,2, ...,m;a = m + 1, ...,n;q = 1,2, ...,n) 

The equations (4.17) together with (4.7) are the required Poincare- 
Hamilton(PH) equations of motion for the nonlinear nonholonomic dynam- 
ical system. They together with the constraint equations (3.3) determine 
the (n + m) quantities x\, X2, x n ; y±, y2, ■■■,y m - In fact, we can find r]j(j = 
1,2, ...,m) as functions of x p , yj and t and then substituting the values of 
r]j into the equations (3.3), determine the remaining r] a (a = m + l,...,n). 
In this way all the ry^'s are determined and then, by using (4.7), all the x p s 
can be obtained as functions of the time t. 



5 Poincare-Cartan Integral Invariant of Nonholo- 
nomic Dynamical System 

We now turn to the investigation of the Poincare-Cartan integral invariants 
for the conservative nonlinear nonholonomic system whose motion is deter- 
mined by the Poincare-Hamilton (PH) equations (4.17) with (4.7). In order 
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to achieve our goal we need to compute, subject to the constraints (3.3) with 
(3.4), the asynchronous variation of the action integral defined by 



S 



f 



Ldt 



(5.1) 



where Ldt expresses the small element of action with L, describing the dy- 
namical behavior of the system, as the Lagrangian function of all the x's, 
r/'s and possibly the time t. 

Let us now consider the motion of the system along the real (actual) 
trajectory between the two positions Pi and P2 corresponding to the initial 
and final instants t\ and t2, respectively. We assume that the varied path is 
determined by means of A-variation in which not only the coordinates but 
also the time varies at the initial and final moments of the motion. Since 
the Lagrangian is a function of class C 2 with respect to all of its arguments, 
using Lemma 2.1, we perform the A-variation of (4.1) according to (2.25) 
to get 

A5= [ 2 (AL + m )dt 
Jti 

which, by virtue of (A. 4) from the APPENDIX, assumes the form 
AS 



,dL 
K dr)j 

rt 2 



J All 



+ 



<2 / 

—f] p Q,o + rj^oX^L - Q L ) dt 



+ 



rti rti 

/ (Lfio)-dt- / 

JU JU 

dL 



'ti 

(Am 



dri P 
d dL* 



ti [dt drjj 



(K*+K k qjVq ) 



dL* 
dvk 



X*L* 



df] a ' 



LOjdt 



(5.2) 



where the superscript " * " shows that the quantities are expressed in terms 
of the independent Poincare parameters of real displacement. 

Since the Poincare equations (4.6) hold along the real(actual) trajectory, 
using (2.11), (2.16) and (4.3), the last integral on the right hand side of (5.2) 
vanishes, and we are left with 



AS 



dL" 
dr]j 



ujj + L*VLq 



\ *2 rti QJ^ q l 

+ I (^—rjpfL + r]^tt X fl L - —rjpfl 

- rjfj,CloXfj,L)dt 



And, in view of (4.10), simplifies to 

AS = (y*io 3 + L*n ) 



11: 



(5.3) 
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Using the relations (2.17) and (4.13), the last result becomes 



AS = {y*nj + H*n )\H (5.4) 

where the Hamiltonian function H* is computed along the actual trajectory 
of the system. Denoting it by — j/q, the asynchronous variation AS of S can 
be put in a more concise form 

AS = yin x \ll (A = 0,l,2,..,m) (5.5) 

If we take G = x x in formula (2.18), we find that 

Ax\ = Vt^X^xx = 0, 1,2, ...,m), (5.6) 

and , by using (2.12) it can readily be shown that are linear combinations 
of the quantities Ax\ in the form 

fi„ = C^Ax A (5.7) 

where ||C^(av)|| is the inverse of the matrix Then the variational 

equation (5.5) becomes 

*2 



AS = (y^AsJ = {y^AxXl (5-8) 

Equation (5.8) yields in concise form the variation of the functional S in the 
space of (2m + 1) variables (x±,X2, ■■■,x m ',y*,y2, ...,y * m ;t) called the "Ex- 
tended Reduced Phase-Space" and may be applied to any functional of this 
type. 

Let us now consider the motion of the system in this space. We note the 
fact that corresponding to the different initial values of the x's and y's we 
obtain a set of initial points at the initial time t\ and through each point we 
can draw the appropriate real paths satisfying (4.17) and giving rise to a set 
of terminal points at the other end of each real path at the same terminal 
time t2- If the initial set of points forms a closed curve C±, correspondingly, 
the terminal set of points will also form a closed curve C2 giving a tube of 
real trajectories. 

Thus integration of (5.5) along a closed curve C (that is, the locus at 
any time t of the points which are initially on C\) that passes from C\ to 
C2, yields 



= £ as = £(ytn x \li) = £(y* x {2) n x w - y *mn x 



ah 
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where J/t , 2/t > ^A^) ^A^ are the values of the y*'s and the f2's at the 
time ti and £2 respectively. In view of the continuity of the motion, the last 
result may be written as 



taken along a closed contour C remain invariant during an arbitrary dis- 
placement with deformation of the system. Thus, we have the following: 

Theorem 5.1. The line integral 



along an arbitrary closed curve C remains invariant with arbitrary defor- 
mation of this curve along the tube of real trajectories of a conservative 
nonlinear nonholonomic dynamical system whose motion is governed by the 
equations (4-17) provided that the relations (3.8) hold. 

We remark that the integral (5.9) is a generalized form of the relative in- 
tegral invariant of Poincare-Cartan given in [3|. It is shown [21 OH] that these 
integrals are important in the study of analytical dynamics of holonomic 
system, since their invariance yields their relationship with the Poincare- 
Hamiltonian system. In what follows, we shall show that this relationship 
can also be established for nonlinear nonholonomic dynamical systems as 
well. 



6 Poincare-Hamiltonian Systems and the Poincare- 
Cartan Integral Invariant. 

In the preceding section we have established the integral invariant of a non- 
linear nonholonomic dynamical system. Now we shall discuss the conse- 
quences of the invariant property of the integral (5.9). Precisely, we prove 




which leads to 





(5.9) 
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the converse of Theorem 5.1 by making use of the property of asynchronous 
variation as discussed in section 2. 

We start with the assumption that the generalized Poincare-Cartan in- 
tegral (5.9) is invariant with respect to the tube of real trajectories of the 
system whose motion is governed by the set of equations of the form 

Vj =^j(x P ,yt,t) (6.1) 

Vj = 4>j{x p ,y* k ,t) (j,k = 1,2, ...,m;p = 1,2, ...,n) (6.2) 

where ipj and <f>j are arbitrary functions to be determined. 

The invariance of the integral (5.9) along the tube of real trajectories, 
in accordance with (2.19), (6.1) and (6.2) implies that 

dl = 

Thus, we have 

= d j>{y*jSlj - H*n ) = j> (dyfoj + yjdSlj - dH*n - H*dn ) (6.3) 

Recalling that relation d(At) = A(dt) it follows that 

H*dfl = H*d{At) = H*A(dt) = A(H*dt) - AH*dt 
Combining this with (2.22), and the identity 

y*Vjdn = y*Vjd(At) = y*Vj^(dt) = A{y*r lj dt) - A{yjrjj)dt, 
we obtain 

-dH*fl - A(H*dt) + AH*dt] 

which takes the form 

= dj> A [(y* jVj - H*)dt] + j> [dy^Slj + y*Ar}jdt - A(y*rjj)dt 

-y*C J qr rj q n r dt - y*C J 0q n q dt - dH*Q + AH*dt 

Since the motion is represented by the closed curve C which is completely 
arbitrary, the first integral on the right hand side vanishes. Thus we have 

= dj> [dypj + y*A Vj dt - A(y* V:j )dt - y]r) q C{ r ^ r dt - y*C J 0q n q dt 

- dH*n + AH*dt] 
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Interchanging the indices r and q and separating the sum over the index 
q = 1,2, n into the sums over the indices i = 1, m and a = m+ 1, n, 
we have 

= d^" [dy*% + y*A Vj dt - A(y* jVj )dt - y* (c 3 riVr + cQ SUdt 

-y*j (c^ ra ri r + cQ n a dt - dH*n + AH*dt 

Using the formulae (2. 18-19) and performing some algebra we get 

= j> [(y* - Ci iV * - C^y*) Slidt + y*A Vj dt - Ay* Vj dt - y*A Vj dt 
-V*j {pla + C™^) (w Q + rj a Q, Q ) dt - rj X a fl dt - rj j X j H*n dt - rj a tt X a H*dt 



3H* BH* 
-Q-rd y *n + fl X H*dt + n 3 X 3 H*dt + Q a X a H*dt + _Ay;dt 



Taking into account (2.16-17) and (3.4), we obtain 



= 



y* - C^y* - Ci iVr y* -y*{Cl a + ^ ra r, r )^ + X t H* + ^X. 11 \ il-.dl 
dH* 



dy* 



dm drji J 

OH* \ 
-Vj(CL + C J ra r]r)vModt - r) a X a H*n dt - — 7 dy*Q + r] a n X a H*dt \ 



dy* 



which, by using (3.11-13), reduces to 



{y* - y*(K> 0i + Kifk) + X*H*} ^dt - Ay*( Vj - —)dt 



(6.4) 



dH* 



As required by the present investigation, we adopt the second viewpoint for 
which 5* f a = for nonlinear nonholonomic dynamical system and multiply 
(3.9) with (£-)* we get 



dL 

drja 



uji = (i = 1, m; a = m + 1, n) 



(6.5) 
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Since the relations (2.17) hold for all the Wj's, we integrate (6.5) along the 
closed curve C and incorporate the result into (6.4), we obtain 

= f [ jtf - y*(Kl t + KM + X*iT - (A?)* (JLj^dt 



^(Vj ~ ^r)dt + \ y*(C«L + Cl aVr )^ m - X*H* Vl 



J 
dL 



dfa 
dr]i 



dH* 



which, in view of the equations (6.1) and (6.2), takes the form 

dL 



= 



- Ay*(^- 



-y*(K> 0i + K> ri r, r )+X*H*-(A?y 
dH* . 



drj 



ay 



dy* 
dL\* 



) + i y*j( C L + CLVr)^Vi ~ X*H*r)„ 



dm 



drja 



dH* 



-y^ + CM^-^y*^ 



dt 



In order that the equations (6.1) and (6.2) must be satisfied for the 
asynchronous variation in which the quantities Jlj and Ay* are arbitrary, 
the coefficient of each of these must vanish. 
This implies that 



-Xfir + (K 3 0i + K 3 ri r h )y* i + (A 



dL 
drja 



; V'i 



dH* 
dy* 



(j = 1,2, ...,m; a = m + l,...,n;r = 1,2, ...,n) 

which yield the desired values of the functions fa and ipi and hence leads to 
the Poincare-Hamilton equations of motion. Also we have 



y]K a + ci aVr )^ Vi - x*H* m + Vl n (A 



dL 
d7] a 



dH* 



Thus the preceding analysis can be summarized in the following: 
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Theorem 6.1. If the line integral (5.9) is invariant under a deformation 

of an arbitrary curve C along the tube of real trajectories, then the motion 

of the conservative nonlinear nonholonomic dynamical system is determined 

by the Poincare-Hamilton(PH) equations (4-17) together with (4 -7), provided 

that the relations (3.8) hold. 

It is remarkable to note that Theorems 5.1 and 6.1 furnish the necessary 
and sufficient condition which allow us to connect the theory of integral 
invariants with the theory of Poincare-Hamiltonian systems. Thus, these 
integrals, as discussed in [2], [3], [121 UHl EHIITT] , are very important in analytical 
dynamics in the sense that they provide another foundation for not only 
holonomic or linear nonholonomic systems but also furnishes a foundation 
for nonlinear nonholonomic dynamical systems whose motion is determined 
by the PH-system of equations (4.17) with (4.7). 



7 A Generalization of The Poincare Linear Inte- 
gral Invariant 

We now turn to a generalization of a theorem analogous to the theorem of 
Poincare [31 [18] as a special case of the results established in the last section. 
To achieve this aim, we proceed with the assumption that the variations 
under consideration are synchronous. This implies that the quantity S7o = 
At = 0. From (2.17) and (2.21), it follows that f2 p = tu p and Ar] p = 8r]p. 
Thus A = d holds and therefore d5 = 5d holds. 

Under such conditions, we may restate Theorem 5.1 as follows 

Theorem 7.1. The line integral 



along any closed curve C that consists of the simultaneous states of the sys- 
tem does not change with arbitrary deformation of this curve along the tube 




(7.1) 
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of real trajectories of the conservative nonlinear nonholonomic dynamical 

system which is described by the PH system of equations (4-17) with (4-V- 

The proof of the theorem adopts a procedure similar to that discussed 
in the preceding sections but without requiring Lemma 2.4. The invariance 
of the integral (7.1) can be obtained by considering the simultaneous states 
of the dynamical system during its motion in the "Reduced Phase-space" of 
2m variables (xj,yj) and using the fact that, contrary to the relation (2.23), 
the d— operation and integration commute with each other. Provided we 
have synchronous variation, we may assert the converse of Theorem 3 in the 
following: 

Theorem 7.2. If the line integral I\, given by (7.1), remains invariant un- 
der an arbitrary deformation along the tube of real trajectories of any closed 
curve C consisting of the simultaneous states of a conservative nonlinear 
nonholonomic dynamical system, then the motion of the system is deter- 
mined by the PH system of equations (4-17) with (4-7). 

We remark that the results above show that the theory of integral invari- 
ants forms another basis for both the Hamiltonian dynamics of holonomic 
systems, linear nonholonomic systems and for nonlinear nonholonomic dy- 
namical systems as well. 

In order to demonstrate the utility of the four theorems, we derive special 
cases of our general results that are analogous to well known results. 

(i) Suppose that all the x's are the Lagrangian coordinates and the r/'s 
are the generalized velocities x's. In this case, the relations (2.17) 
reduce to the result given in [25] and [26] by Vujanovic and Whit- 
taker, respectively; the operators Xq and X p 's becomes and ^-'s. 
Consequently all the Cq 's and C^'s vanish. Theorem 7.1 and 7.2 
furnish the results analogous to those obtained in [7], [12], [18], and 
|26j by Cartan, Gantmacher, Pars and Whittaker, respectively while 
Theorems 7.1 and 7.2 subsume the results that are discussed in |19j 
and [20] by Poincare and in [18] by Pars. It is to be noted that our 
results are analogous to these but the content is quite different, since 
our system is nonlinear nonholonomic. 
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(ii) If the group variables are the quasi-variables (nonholonomic coordi- 
nates) 7r's then the 77's becomes 7r's and the relations (2.10) express 
these as non-integrable linear combinations of the quasi- velocities and 
all the Cg p 's reduce to Hamel-Boltzmann's three indexed symbols 7^ . 
In this case our theorems subsume the results analogous to those ob- 
tained in [9] by Djukic. 



APPENDIX 

Asynchronous Variation of Functional (5.1) 

Performing the A-variation of (5.1) according to the Lemma 2.4, we have 

eta 



AS = [ 2 (AL + Ln )dt 
Jti 



ft2 

= / (— A Vp + %X^L + (Lfi ) - Ln )dt 
Jti orjp 

where we have used Def. 2.3 and the identity (LClo) = LCIq + L£Iq. Taking 
into account (2.16), (2.17) and (2.21), the variation AS of S becomes 

f t2 dL f t2 dL 

^ s = ( o— 8Vp + oj p X p L)dt + I (—rjpQo + "H^IqX^L - VL L)dt 
Jti orjp J tl orj p 

+ [\m )dt (A-l) 
Jti 

Let 

<' t2 , dL 



f 2 dL 

1=1 {-^—$Vp + u p XpL)dt, 
Jti "Vp 



separating the sum over the index p from 1 to n into the sums over j from 
1 to m and over a from (m + 1) to n and using the equation of motion (4.1) 
to get 



t2 



dL* dL d<b a \ , dL 



h L v drjj dr] a drjj J drj, 



Srjj + — — Srj a + ujpXpL \ dt. 



Here the '*' over the quantities show that they are expressed in terms of the 
independent parameters r\j of real displacement. 
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Following the second viewpoint and using relation (2.20), we have 



* 2 dL* 



i dr lj 



(uj + C 3 0q u q + C 3 qr r] q u r )dt 
rt * f dL 



* 2 dL 



h dri a drjj 



+ C 3 qr r] q u r )dt + J | — (cj Q + Co q uj q + Cg r r] q uj r ) + j dt. 

Integrating by parts, the first term of each integral on the right hand side 
of the last result, we get 



dL* 



dL 



t2 d .dL*. 



r 2 dL* 



drj a drj - ' 



* 2 Z** 2 d , dL 
+ to 

tl Jtl 



3 dt drj a drjj 



. dL , 
dr] a 



*2 /-t 2 



d , dL 



tl Jtl dt 9r? a Jtl 



t 2 



(C > 9 + C^uv) + w p X p L dt 



which by use of (3.4), reduces to 

"* 2 f d ,0L* . 



Or/-,- 5 
<9L d 



) -r^< 



<9L* 



d , dL d(f) a dL 



<^(— )^-) + 



dt x d-q a ' drjj ' ' Or/a <9?7j 



drj a dt drjj 

Again breaking the sum over the indices q and r from 1 to n into the 
sums over the indices j, k from 1 torn and over the indices a, (5 from (m+ 1) 
to n, the last expression takes the form 



I = 



dL* 
dr] : 



dL d 



'I] "'13 

dL d<p a 



+ c^ a ) - *3 WaJt ^) + ^ ^-(^ + cy*t + c 3 qkV , 

dL 1 



Interchanging the indices j and k in the second and fourth terms of the 
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integrand, using (3.4) and rearranging the terms, we find that 



dL* 
« — w i 



1 2 



tl J tl 



' — ( — ) — 



{C 0j + C O a -Q^) + Vq(C qJ + C qa-Q^) 



dL d d<f> a dL dcf) a 
3 dr] a dt drjj dr) a dr] k 



( c oi + c k w-E-) + + c W~^)m \ 



d(j>8. 



dq 



dL 



dr\ a 



d<j>0 



drij 



dt 



Taking into account equations (3.11) and (3.12), the last result reduces to 



/ = 



dL* 



*2 



Ju 



d ,dL 



dL 



H 3 

(—Y f _ ( d( ^ a \ 

Wi %J \dt ( d Vj ) drj k 



CO.; 



h Jt! v dt d m dm 



(A-2) 



— UpXpLj dt 

where the "*" over the quantities (J^-)* indicates that they are expressed 
in terms of the independent parameters rjj. 

Let us consider the term u> p X p L which, in view of (4.4), can be written 

as 

dL 

ojpXpL = ujpXpL — Up— — X p (j) a 

orjot 

= LOjXjL* + LO a X a L* — LOj— Xj(f) a — LOp- Xp(j) a 

OTj a OT] a 

9L V A, 
L)j- Aplpa 



* d<p a dL 

drjj drj a J drjj J dr], 



where to obtain this result we have separated the sum over the index p = 
1,2, ...,n into the sums over j = 1,2, ...,m and a = m + 1, ...,n and also 
used the relation (4.3). Simplifying the last expression, we get 



ujpXpL = ujjiXj + -j^j.\„)i: - ujj — iXj + '-jj-j.X ,)<■>„. 
which together with (3.13), becomes 

dL 

ujpXpL = UjX*L* - — )*X*4> a , 

OT] a 



(A-3) 
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where we have expressed all the quantities in terms of the independent pa- 
rameters rj/s of real displacement. This allows us to write (A.l) as 



dL* 
dL 



d .d<t> 



dt drjj drjk 



+ K k qjVq ) + {K& + K" Vq ) - X* 



ujjdt 



dL* 
Or), 



-co, 



tl 



tl Jti V dt d n. 



d ,dL* 



Pit* 

-(^+K k qjVq ) Wk -X*L*-(A ? ) 



dL 



dri a 
(A-4) 



which expresses the asynchronous variation of the action integral (5.1). 



ujjdt 
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